On Particle Transport and Radiation Production in Sub-Larmor-Scale Electromagnetic Turbulence 
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The relation of particle transport of relativistic particles in plasmas with high-amplitude isotropic sub-Larmor- 
scale magnetic turbulence to the spectra of radiation simultaneously produced by these particles is investigated 
both analytically and numerically. We have found that in the asymptotic regime of very small particle deflections 
the pitch angle diffusion coefficient is directly related to the spectrum of the emitted radiation. Moreover, this 
spectrum provides much information about the statistical properties of the underlying magnetic turbulence. The 
transition from small- to large-scale jitter to synchrotron radiation regimes as a function of turbulence properties 
has also been explored. These results can readily be used to diagnose laboratory and astrophysical plasmas. 



I. INTRODUCTION 

High-amplitude sub-Larmor-scale electromagnetic turbu- 
lence is ubiquitous in high-energy density environments, such 
as laboratory plasmas produced by high-intensity lasers, e.g.. 
National Ignition Facility, Omega, Hercules, Trident, and 
others fT^HI, and in astrophysical and space plasmas, e.g., 
at high-Mach-number collisionless shocks in weakly magne- 
tized plasmas lEJd, upstream regions of quasi-parallel shocks 
jB., 9|, sites of magnetic reconnection IfTOl [TT] and others. 
Studies of plasmas and turbulence in these environments are 
important for fusion energy sciences and the inertial confine- 
ment concept [1 , 4|, in particular, as well as to numerous as- 
trophysical systems such as gamma-ray bursts 1. 12.-. 14 J . su- 
pernovae blast waves ifTSl . jets of quasars and active galactic 
nuclei |161, shocks in the interplanetary medium fTTl, solar 
flares 1 18 1 and many more. Such turbulence can be of various 
origin and thus have rather different properties, from being 
purely magnetic (Weibel) turbulence fT9','20l, to various types 
of electromagnetic turbulence (for example, whistler wave tur- 
bulence or turbulence produced by filamentation/mixed mode 
instability), to purely electrostatic Langmuir turbulence 1,211 - 

Despite substantial differences, these turbulences share one 
thing in common: the electromagnetic fields vary on scales 
much smaller than the characteristic curvature radius of the 
particles' paths (i.e., the Larmor radius in most cases), which 
is typically of the order of the plasma inertial length (skin 
depth), so that the particle orbits are never a well-defined Lar- 
mor circle. If the electromagnetic fields are random, which is 
usually the case of turbulence because of the random phases of 
fluctuations, their paths diffusively diverge due to pitch-angle 
diffusion. Radiation simultaneously produced by these parti- 
cles is neither cyclotron nor synchrotron (for non-relativistic 
or relativistic particles, respectively) but, instead, carries in- 
formation about the spectrum of turbulent fluctuations. In this 



paper we explore the relation of the transport of relativistic 
particles in isotropic three-dimensional magnetic turbulence 
and the radiation spectra simultaneously produced by these 
particles. 

We have found that in the small-angle-deflection regime, 
the radiation spectrum agrees with the jitter radiation predic- 
tion lfT3l 124^ 271. We also demonstrate that the pitch-angle 
diffusion coefficient is directly related to and can readily be 
deduced from spectra of the emitted radiation. This provides 
a unique way to remotely diagnose high-energy-density plas- 
mas, both in laboratory experiments and in astrophysical sys- 
tems. This can be particular interest, for example, for the 
physics of collisionless shocks and Fermi particle accelera- 
tion in them, because the latter is intimately related to the dif- 
fusion coefficient in the upstream and downstream regions of 
such shocks. 

The rest of the paper is organized as follows. Section [II] 
presents the analytic theory. Sections III and IV describe the 
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numerical techniques employed and the obtained simulation 
results. SectionlVlis the conclusions. 



n. ANALYTIC THEORY 

Let us consider a relativistic electron moving through a 
non-uniform, inhomogeneous, random, small-scale magnetic 
field (and we assume that this magnetic micro-turbulence is 
statistically homogeneous and isotropic). Because of the ran- 
dom Lorentz force on the electron, it's acceleration and veloc- 
ity vectors vary stochastically, leading to a random (diffusive) 
trajectory. We define the field turbulence to be "small-scale" 
when the electron's Larmor radius, pe — •jPnieC^ /eB± 
(where 7 is the electron's Lorentz factor, /3 = v/c is the di- 
mensionless particle velocity, nie is the electron mass, c is the 
speed of light, e is the electric charge, and B± is the compo- 
nent of the magnetic field perpendicular to the electron's ve- 
locity vector) is greater than, or comparable to, the character- 
istic correlation scale of the magnetic field, Xb, i.e., Pe <^ As. 
Hereafter, we consider ultra-relativistic particles only, j ^ 1. 

For small deflections, the deflection angle of the velocity 
(with respect to the particle's initial direction of motion) is ap- 
proximately the ratio of the change in the electron's transverse 



momentum to its initial momentum. The former is ~ Flt\, 
where Fl = (e/c) v x B is the transverse Lorentz force, and 
Tx is the transit time, which is the time required to traverse the 
the scale of the field's inhomogeneity, i.e., the field correlation 
length, Xb- For an ultra-relativistic particle (whose Lorentz 
factor is 7 ;:^ 1), t\ ^ Xb/c. The particle's total momentum 
is, likewise, p ^ jm^c. The change in the transverse momen- 
tum is thus, p± ^ Flt\ ^ cB^Xb/c. An ultra-relativistic 
electron will only experience small deviations to its original 
path; consequently, the deflection angle over the field corre- 
lation length will be a\ « P^_/p ~ eB^XB/^irn^c^. The 
subsequent deflection will be in a random direction, because 
the field is uncorrected over the scales greater than As, hence 
the particle motion is diffusive. As for any diffusive process, 
the mean squared pitch angle grows linearly with time as 

{a") = D^at. (1) 

The pitch-angle diffusion coefficient is, by definition, the ratio 
of the square of the deflection angle in a coherent patch to the 
transit time over this patch, that is 



Ae ~ l/y 
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where an average square magnetic field, (i?^) has been sub- 
stituted for B^. Note that the diffusion coefficient depends 
on both statistical properties of the magnetic field, namely its 
strength and the typical correlation scale. 

Next, the diffusing particle experiences accelerations and, 
hence, produces radiation. The radiation from an ultra- 
relativistic charge is beamed in a cone with a narrow angle 
of A6 ^ 1/7. The ratio of the particle deflection angle to the 
beaming angle is Il24ll25l : 



a 



cB^Xb 



2tt :-— = 6. 

meC^KB 



(3) 



where we used that B±^ ^ (B^)^'^, as before and the charac- 
teristic wavenumber of the turbulence, i.e., k at which |BkP 
has the maximum, is fc^ ^ 27rA^^. Note however that the 
jitter parameter is scale-dependent, therefore at each scale k, 
the jitter parameter is 
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where (Bf)^^"^ is the rms field at scale k, i.e., it is average 
over a range of scales from fc to fc + Ak with A/c ^ k, that is 
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The radiation spectrum and diffusion may, generally, be sen- 
sitive to this fc-local value Sk, though the overall spectrum and 
diffusion are indeed determined by its averaged value, S. One 
important conclusion can readily be made from comparison 
of Eqs. (J2]) and ([3]l. The jitter parameter and the pitch-angle 
diffusion coefficient are directly related as 
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Figure 1. Radiation regimes, (a) Large-angle jitter regime, 1 < 5 < 
7; radiation is only seen along certain segments of the particle path, 
thus resulting in the spectrum that is synchotron-like at and above the 
peak but differing from synchrotron at low frequencies, (b) Small- 
angle jitter regime, 5 < 1; radiation is seen from the entire trajectory, 
thus the spectrum depends on the underlying spectrum of electro- 
magnetic turbulence. 



There are several regimes of interest. First, (5 — > 00 regime 
corresponds to the classical synchrotron radiation regime: the 
particle orbit is a circle in the plane orthogonal to the mag- 
netic field and the field itself is homogeneous. Second, the 
regime with (5 > 7 is very similar to synchrotron, but the par- 
ticle's guiding center is slowly drifting in an inhomogeneous 
magnetic field. The produced spectrum is well represented by 
the synchrotron spectrum in a weakly inhomogeneous mag- 
netic field and it slowly evolves in time due to particle diffu- 
sion through regions with different field strength. This regime 
can be referred to as the diffusive synchrotron regime. Third, 
when 1 < (5 < 7, the particle does not complete it's Larmor 
orbit because the B-field varies on a shorter scale. In this case, 
an onlooking observer would see radiation from only short in- 
tervals of the particle's trajectory (i.e., whenever the trajectory 
is near the line-of-sight), as in synchrotron, but these intervals 
are randomly distributed. This is the case of the large-angle 
jitter regime. The radiation is similar to synchrotron radia- 
tion near the spectral peak and above, but differ significantly 
from it at lower frequencies, see Ref. [25| for details. Fi- 
nally, if (5 < 1, the particle's deviations are extremely small 
compared to its beaming angle, so the radiation seen by an ob- 
server is produced over a large number of incoherent patches 
of the i?-field, i.e., almost throughout the entire trajectory of 
the particle III31 1241 125 L as is shown in Fig. [T]. Thus, the re- 
sulting radiation markedly differs from synchrotron radiation, 
although the total radiated power of radiation, Ptot = dW/dt, 
produced by a paiticle in all these regimes, e.g., jitter and syn- 
chrotron, are identical: 
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where r^ = e^ /meC^ is the classical electron radius and tUc ~ 
^'^e{B^)^^'^ /irieC is the characteristic synchrotron frequency 
of radiation in a uniform or weakly inhomogeneous (on the 
Larmor scale) magnetic field. 

Radiation produced in the small-angle jitter regime, S < 1, 
is the most sensitive to and provides the most information 
on the turbulence properties, hence, we mostly consider this 
regime in our paper. The characteristic frequency of radiation 
can be estimated by considering the virtual photon approxi- 



mation. The correlation scale A^, as seen in the electron's 
rest frame, is related to Xb in the lab frame by the Lorentz 
transformation (i.e., A^ = Xg/'y). The electron's transverse 
acceleration occurs over the length scale of Xb, which corre- 
sponds to a characteristic time scale, in the electron's frame, 
of r' ^ '^'s/c. Thus, the radiation emitted by the electron 
has a characteristic frequency, in the electron's rest frame, 
of uj' ^ 2tt/t' ^ 2TTc/X'g ^ 2tt^c/Xb- Transforming 
back to the lab frame picks up another factor of 7 from the 
Lorentz transformation, giving the characteristic frequency of 



27rc7 /As, or: 



LOj ^ 7 kBC ~ inuJc/S. 



(8) 



The spectral power at the characteristic (spectral break) fre- 
quency is 
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Note that this equation allows one to determine the jitter pa- 
rameter directly from spectral observations: 
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From Eqs. (J6]l, (J8]l one obtains the relation between the 
spectral break of radiation and the pitch-angle diffusion coef- 
ficient: 
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One can also relate the diffusion coefficient to the spectral 
power at the break frequency from Eq. (J9|l. We obtain 
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We emphasize that Daa is, thus, a directly measurable quan- 
tity because it solely depends on the radiation spectral power 
at the break frequency, P{ujj), and the Lorenz factor, 7, of 
the radiation-emitting electrons. Thus, radiative diagnostics 
should be a very useful technique to study laboratory plasmas 
and explore plasma conditions in astrophysical plasmas. 

In general, turbulence has a range of scales rather than a sin- 
gle scale kB- We explore how spectral and transport properties 
are modified in this more realistic case. In the analysis, we 
assume the isotropic three-dimensional magnetic turbulence 
with a power law turbulent spectrum: 
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if fcmin < fc < fcmax and zcro otherwise. Here Bk is the spatial 
Fourier transform of the magnetic field, k is the wave vector, 
and index /i is a positive number. Notice that we consider a 
static isotropic turbulence, hence the field distribution is inde- 
pendent of time and the wave vector direction. Note that such 
magnetic turbulence is a natural outcome of the non-linear 
Weibel/filamentation instability, which occurs at relativistic 



collisionless shocks and in laser-produced plasmas. It has 
been shown lfT3l |25| - |27J that monoenergetic relativistic elec- 
trons in such turbulence produce flat angle-averaged spectra 
below the spectral (jitter) break and power-law spectra above 
the break, that is 

if LO < LUj , 

if ujj < UJ < ujb, (14) 

if ujb < UJ, 

where the spectral (jitter) break is defined in this case as 




YK 



(15) 



because the field Bk at k^nm has the most power (for /i > 0) 
and, hence, kj„i„ plays the role of fc^ in Eq. ([8]). Similarly, the 
high-frequency break is 



uJb = 7^M™xC. 



(16) 



Thus, by observing the radiation spectrum, one can probe the 
structure of the magnetic turbulence that generates it. Note 
that the angle-averaged single electron spectra are equivalent 
to the ensemble-averaged spectra for isotropic particle distri- 
bution functions. We will use the latter in our numerical study. 



III. NUMERICAL MODEL 

Our goal is to explore, via simulations of particle dynamics 
in magnetic turbulence, the diffusive and radiative properties 
of plasmas and how they are related. We perform our sim- 
ulations from first-principles. We use a 3D simulation box 
with periodic boundary conditions in all directions. Relativis- 
tic electrons are test particles moving in the preset magnetic 
fields but do not interact with each other, nor do they induce 
any fields. Radiative energy losses are considered negligible 
compared to the energies of individual particles, hence ne- 
glected. Motion of each electron is, thus, solely determined 
by the Lorentz force equation given by: 
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(17) 



where (3 = \-{t)/c. Note that this Lorentz force conserves 
particle's energy, therefore the magnitude of the velocity is 
constant in time and only the direction is changing. The sim- 
ulation can be divided into two principle stages (see Ref 1281 
for a detailed description of the numerical implementation). 
First, the turbulent magnetic field is generated from a given 
spectral distribution in Fourier space. This field is created on 
a lattice that is then interpolated, so that a "continuous" field 
is represented. The interpolation was implemented by way of 
divergenceless matrix-valued radial basis functions (see Ref. 
129], for a discussion). This method begins with a radial func- 
tion - in our case, one of the simplest, 0(r) = e^'^^ (where 



e is a scaling factor, and r 



c- + y + z ). Then, a set 



of divergence-free matrix-valued radial basis functions is ob- 
tained from the transformation |29|: 



$(r) = (VV^-l3x3V2)0(r), 



(18) 



where VV^ is the second-order, 3 x 3-matrix differential op- 
erator and I3x3 is the 3x3 identity matrix. These interpolants 
were applied to the interior of each lattice "cube" (i.e. the 
space between a single lattice point and the five immediately 
adjacent points). The second stage in our model then involves 
the numerical solution of the equation of motion for each par- 
ticle, from which {a^) and the radiation spectra are obtained. 
We will first turn our attention to the generation of the mag- 
netic field. 

Generation of the magnetic field distribution is more con- 
venient in Fourier space for two reasons. Firstly, it is simpler 
to specify a particular spectral distribution in Fourier space 
directly, rather than attempting to find a field in real space 
that would fit the chosen spectral distribution. Secondly, any 
physically realizable field should satisfy Maxwell's equations, 
thus its divergence must be zero. Producing a divergenceless 
field in Fourier space is far simpler than generating one in real 
space. This is because Gauss' law, V • B = 0, is an algebraic 
equation in Fourier space, k-Bk — 0, for each Fourier compo- 
nent, Bk, corresponding to the wave vector k. Although our 
code can handle a wide variety of magnetic spectral distribu- 
tions, this study is limited to the case of isotropic magnetic 
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turbulence, described in Eq. ( 13 1, and we have left considera- 
tion of more sophisticated models for the future. 

The next stage in the simulation starts with the numerical 
solution of the equation of motion, Eq. (fTTl). This was done 
via a fixed step 4* -order Runge-Kutta-Nystrom method. With 
all the particle positions, velocities, and accelerations calcu- 
lated, the radiation spectrum is obtained from the equation 
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(19) 

which represents the radiative spectral energy, dW per unit 
frequency, duj, and per unit solid angle, dil. In this equation 
r{t) is the particle's position at the retarded time t, n is the 
unit vector pointing from r(t) to the observer and /3 = df3/dt. 
There is some freedom in the calculation of the total radiation 
spectrum. One can add the spectra coherently (i.e., by taking 
the Fourier transform of the "summed over" radiation fields 
of each particle). In this case, only a single integration would 
be needed. Alternatively, we can add the spectra incoherently 
(i.e., by integrating each particle's radiation field separately, 
and then summing the results of each integration). Both ap- 
proaches result in the same spectra, but those obtained in the 
former approach are noisier for the same number of simulation 
particles. Hence we use the second approach in our study. 

Our code was tested in various set-ups. In one, we con- 
sidered motion of a single relativistic particle in a uniform 
magnetic field. Figure l2] shows the obtained radiation spec- 
trum and the corresponding theoretical synchrotron spectrum. 
Here, the blue curve represents the numerically resolved syn- 
chrotron harmonics (which are integer multiples of the gy- 
rofrequency, ujb = eB±/meC). We see excellent agreement 
with the analytical result, indicated in red. 

We also verified that motion of a particle in small-scale ran- 
dom magnetic fields results in pitch-angle diffusion and the 
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Figure 2. Radiation spectrum (i.e., the total radiated energy, dW 
per unit solid angle, dil per unit frequency, duj vs. frequency, uj) 
of a single relativistic charge moving through a uniform magnetic 
field. The numerical solution is indicated in blue, the red line is 
the analytical solution, and they agree very well. The spectrum is 
peaked at the synchrotron frequency ojc = (3/2)7^a;_B, where cjb = 
eB±/'fmsC is the electron gyrofrequency. In this and other spectral 
plots, the radiation power is arbitrarily normalized. 
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Figure 3. The trajectories of 50 monoenergetic (7 = 3) particles 
through a turbulent magnetic field (5 ~ 1) projected on to the x- 
y plane. Each particle (denoted by a unique color) starts from an 
origin with a random initial velocity. The axes are x and y positions 
in simulation units. 



paiticle orbit is chaotic. The chaotic nature inherent in the 
particle motion is illustrated in Figure [3] Here, 50 monoener- 
getic (7 = 3) particles are sent out from an origin in random 
directions. The variabiUty in the particle motion is seen after 
shortly leaving the origin. The diffusive nature of the particle 
motion is evident from Figure [4] and Figure [5] In Figure |4] 
the velocity space of a single relativistic particle (7 = 5) is 
plotted over a very large simulation time. Notice that the par- 
ticle velocity is confined to a sphere (i.e., 7 is constant), and 
that the velocity vector diffusively visits various directions in 
the course of the particle's motion. In Figure l5] the average 
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Figure 6. Pitch-angle diffusion coefficient, Daa vs. particle inverse- 
square Lorentz factor, 1/7^. The "blue squares" indicate the Daa 
obtained directly from simulation (as the slope of (a^) vs. time), 
while the "red triangles" are the analytical, given by Eq. \\\\ , pitch- 
angle diffusion coefficients. 



Figure 4. Velocity space of a single particle (7 = 5) moving through 
isotropic magnetic turbulence {5 ~ 1). The axes are the components 
of the velocity, which are in units of c. Notice that, although the ve- 
locity vector of the particle varies randomly (and, over enough time, 
visits all possible directions), its magnitude is constant. 
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Figure 5. Average square pitch-angle vs. time (in simulation units). 
Here, Daa is the numerically obtained pitch-angle diffusion coeffi- 
cient (i.e., the slope of the line). The slopes are estimated via lines of 
best fit (indicated by black dashed lines). Notice that, although these 
results differ by orders of magnitude in the slope (and, additionally, 
particle number), there is not an appreciable deviation from linear 
behavior (as long as 5 <C 7). 



square pitch-angle, {o?) is plotted versus i for three different 
cases. The linear dependence is indicative of pitch-angle dif- 
fusion, cf. Eq. ([1]). In all three cases, (5 w 1. Note that the 
diffusion approximation of these particles' motion is accurate 
so long as (5 <C 7. 

Finally, we verified the numerical convergence of the re- 



sults by changing the size of the simulation box and the to- 
tal number of simulation particles. Because of vastly differ- 
ent computational requirements in various cases, however, we 
were sometimes forced to use a smaller number of particles 
than in other scientific runs. This introduced more noise in 
some spectra presented below. 



IV. NUMERICAL RESULTS 

In Section [n] we made a number of theoretical predictions 
concerning radiation and transport properties of plasmas with 
small-scale turbulent magnetic fields. Here we check these re- 
sults with numerical simulations and then farther explore how 
radiation spectra depend on underlying magnetic field distri- 
butions. 

First of all, we explore how the pitch-angle diffusion coef- 
ficient depends on various parameters, cf. Eq. Q, namely the 
particle's Lorentz factor, 7, the magnetic field strength (B^) 
and the jitter parameter 8. 

Figurel6]shows the 7 dependence oiDaa- Close agreement 
between the theoretical and numerical diffusion coefficients 
with varying Lorentz factor is evident. Next, D^a depends on 
two additional parameters, namely /;;,„;„ and (S^). Instead of 
plotting Daa vs. /cmin directly, we chose to plot the diffusion 
coefficient vs. the frequency (Figurel7| of Wj - which is given 
in Eq. (|8]l - in the radiation spectrum. 

Figure IT] represents cases with different values of /cmin's, 
fcmax's, Ai's, A'p's, {B'^)'s, and 7's, which is why the points 
show large spread (the parameters used for each data point are 
listed in Table ll|. The dependence of Daa on the magnetic 
field strength, (B^), is illustrated in Figures] which agrees 
with the prediction given by Eq. (J2]). It is worth noting that 
the apparent difference between the theoretical and numerical 
diffusion coefficients in these plots is due to the approximate 
nature of our theoretical analysis in Section [III 
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Figure 7. Daa vs. the frequency of Uj, i.e., the frequency at which 
the radiation spectrum peaks: Eq. Jsl (in this case, obtained numer- 
ically). Once more, the "blue squares" indicate the Daa obtained 
directly from simulation while the "red triangles" are the analytical 
Dec, given by Eq. ijTTJ. This figure represents cases with various 
values of fc,^i„, fc„,ax, Np, (B^), Ai, and 7, hence the large spread 
in Daa& (see Table [l] for a listing of the parameters). Notice that, 
although there is considerable variation between the difference in 
the analytical and numerical diffusion coefficients, the difference is 
never greater than a factor of a few. 



10- 



10' 



10' 



10' 



Np = 10,000 particles 
analytical: triangles 
numerical: squares 



lo-'n 

10^ 



10' 



lo-* 

<B^> 



10- 



10" 



Figure 8. Daa vs. {B^). As before, the "blue squares" indicate the 
Daa obtained directly from simulation while the "red triangles" are 
the analytical Daa, given by Eq. \\l). 



We now present radiation spectra and demonstrate how they 
are related to the parameters of the turbulent magnetic field. In 
all runs, we specify the field strength with the jitter parame- 
ter, given by Eq. ([3]). Also, we will plot the angle-integrated 
spectra dW/du) rather than dW/dVtdtu. The former represent 
the spectra from an ensemble of particles with an isotropic 
velocity distribution. Thus, by summing the individual spec- 
tra of each particle, the solid angle, dO, has been effectively 
"integrated" out. 

Figure l9j shows the radiation spectra as a function of the 



#5 Ai 7 /i fc„ 



{5^)1/2 Np 



1 0.12 0.0025 5 -3 1.3 32.2 0.024 4000 

2 0.63 0.0100 8 -3 1.0 16.1 0.100 2000 

3 0.47 0.0100 7 -3 0.6 16.1 0.047 500 

4 0.47 0.0100 3 -3 0.6 16.1 0.047 500 

5 0.94 0.0100 5 -3 0.3 16.1 0.047 500 

Table 1. Table of parameters used in Figure [7] 
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Figure 9. Radiation spectrum (dW/duj vs uj) illustrating the fcmin 
dependence. Recall that fc,„i„ is the characteristic wave number of 
the turbulent magnetic field (given by a negatively sloped power law 
spectral distribution). The scaling with fcn,i„ agrees with Eq. iflSj. 



normalized frequency, lo/loj for two different values of fcmin- 
Other parameters are: 7 = 5, /i = 3, (i?^) = 0.047, and 
fcmax = 1287r/25 (in both runs), 5 = 0.47, iVp = 500 (in 
the fcmin = 7r/5 case), 5 = 0.94, and Np = 2000 (in the 
fcrain = tt/IO casc). The spectral break occurs somewhat be- 
low LUj but it correctly scales with fcmin, cf Eq. (15i. The 



spectrum is flat below the spectral break and falls-off above it 
roughly as lo^^, which is expected for the used value of /i = 3. 
There is a hint of the second break at the high-frequency end 
due to the high-A: cut-off of the magnetic spectrum. To demon- 
strate this we performed two long runs with different values 
of /cmax/fcmin = 25.6, 51.2 (othcr parameters being the same). 
Figure [T0| shows the results of these runs. One can see that 
the extent of the power-law part of spectrum increases with 
increasing value of fcmax/fcmin- The frequency in the figure 
is normalized by fc^ax rather than fcmin and the alignment of 
the high-frequency breaks is evident. Thus, this break is in- 
deed due to the small-scale cut-off of the magnetic field spec- 
trum, cf Eq. (fT6]l. The spectral slope in the power-law region 
between the two breaks is determined by the magnetic field 
spectral index ^. Figure[TT]shows the results of two runs with 
two values of /i = 2.5, 3. Other parameters, for both cases, 
are: 7 = 5, fc„,i„ = 7r/5, k„iax = 2567r/25, Np — 4000, and 
(S^) = 0.047. One sees that the spectral power law index is 
indeed close to the value of —^ + 2, as is given by Eq. ( [T4| . 
Thus, the radiation spectrum in the small-angle jitter regime 
is largely determined by the B-field spectrum. 
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Figure 10. Radiation spectrum illustrating tlie femax dependence (the 
two spectra differ in fcmax by a factor of 2). The transition from a 
power law to a steep drop off occurs at ojt ~ 7^fcmaxC. 
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Figure 1 1 . An illustration of the magnetic spectral index, fi depen- 
dence. Given a magnetic spectral distribution |-Bfep = k~'^, where 
/i is positive, the expected slope (in a log-log plot) of the particle ra- 
diation spectrum is — /i + 2; following ljj = 7^fcminC, it will be zero 
for lesser frequencies, in agreement with Eq. ( |14^ . 



The spectrum also depends on the particles' energy. Figure 
10 demonstrates that the break frequency Wf, scales ex 7^ and, 
moreover, the shape of the radiation spectrum is independent 



of the particle's Lorentz factor, as seen in Figure 12 



Finally, the radiation spectrum depends on the jitter param- 
eter, S and on the regime of radiation, small-angle vs. large- 
angle, i.e., (5 < 1 vs. (5 > 1, respectively. In Figure 13 
we have superimposed spectra of increasing jitter parameter. 
We see that radiation is in the small-angle jitter regime for 
i5 < 1, which is evident from the overall spectral shape: a 
flat part below the main break, a decaying power-law above 
it, and independence of the spectral slopes and the positions 
of the break on the value of (5-parameter. As 6 increases to 
values somewhat greater than unity, 1 < (5 < 7, the spectrum 
acquires some synchrotron features: the characteristic slope 
of 1/3 right after the peak and the change of the peak posi- 



Figure 12. Spectra for 7 = 5 and 7 = 8 superimposed. By nor- 
malization of the frequency by Uj = 7'^fcminC, one can clearly see 
the invariance of the spectral shape with respect to 7. Additionally, 
there is a striking contrast in the degree of "noise" in the two spectra. 
The two spectra differ in particle number by a factor of 2, with the 
sharpest spectrum given by the 7 = 5, with Np = 8000. 
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Figure 13. Radiation spectra with variable jitter parameters. As can 
be seen, if the jitter parameter, 5 is less than one - the radiation 
spectrum is restricted to the small-angle jitter regime. As 5 becomes 
greater than one, a transition occurs to the large-angle jitter regime, 
and then the synchrotron spectrum (with its distinctive 1/3 slope) 
emerges as 5 becomes much greater than one. The Lorentz factor is 
5 in each case. 



tion with 6 or, equivalently, with (-B^), indicating that it is at 
the synchrotron frequency ujc- However, the high-frequency 
power-law is still well-established above the spectral peak - 
in contrast to the synchrotron exponential decay at large w's. 



The jitter break at ojj 



: (see Ref. ["ZSl) is not well seen 



because of a rather narrow "dynamical range" of the spectrum 
between ojc and the fundamental frequency ^ oJc/j^ for our 
7 =: 5 runs. At even higher values of the jitter parameter, 
(5 > 7, the spectrum becomes purely synchrotron. 

The regime of moderately large deflections, 5 > 1, is inter- 
esting, as it produces a "hybrid" jitter-synchrotron spectrum. 
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Figure 14. Spectrum from the small-scale field contribution (5k G 
[0.1, 1]) superimposed with the spectrum from the large-scale field 
contribution (S^ £ [lj4]), their sum, and the full scale interval 
(Sk £ [0.1, 4]). Each wave number of the magnetic field gives rise 
to its own "scale", and therefore its own jitter parameter. The small- 
scale field contributes to a jitter spectrum (in red), while the large- 
scale field contributes to a synchotron-like spectrum (in blue). The 
orange dashed line indicates a typical synchrotron spectral shape. 
The fact that the sum of the spectra from the two scales (purple 
curve) entirely agrees with the spectrum from the full 5 interval 
(green curve) illustrates the independence of radiation production at 
different scales. 



In this regime, S^ of Eq. (J4| exceeds unity for some k and 
below unity for others. In this regime, the global (5-parameter 
gives a poor description of the radiation regime, since a low-fc 
part of the magnetic field power spectrum is in the large-angle 
regime, whereas a high-fc end produces small-angle jitter radi- 
ation. To demonstrate this, we have done the following experi- 



ment. In Figure 14 we have divided a wave number range into 
the small-scale portion (i.e., all k for 6k < 1) and the large- 
scale portion (i.e., all k for 6k > 1). We see that the large-scale 
component produces a synchrotron-like spectrum (blue, solid 
curve); the dashed orange curve designates the typical syn- 
chrotron spectral shape, and the small-scale portion is in the 
small-angle jitter regime (red, solid curve). Additionally, we 
see that the sum of the two spectra (purple, dot-dashed curve) 
is nearly identical to the spectrum obtained from the entire 
wave number interval (green curve). 



V. CONCLUSIONS 



angle diffusion coefficient and the simultaneously produced 
radiation spectrum are tightly related, cf. Eqs. Q and (|6|. 
Moreover, the diffusion coefficient and the jitter parameter of 
turbulence, 6 = Aa/(l/7), can readily be determined from 
the spectral information alone, cf. Eqs. (Ill, ( 12 1 and ( 10 1, 



respectively. These theoretical results have further been con- 
firmed with first-principles numerical simulations. 

Furthermore, we numerically explored how much informa- 
tion about the statistical properties of the underlying magnetic 
turbulence is present in the radiation spectrum produced by 
relativistic particles propagating though it. We have found 
that in the small-deflection-angle regime the relation of the 
radiation spectrum to the Bk power spectrum is very tight. 
For example, for the power-law spectrum of i3-field, Eq. 
( 13 1, the produced radiation spectrum is described by a dou- 
ble power-law (even for monoenergetic electrons) — in clear 
contrast with the standard synchrotron theory. The low- and 
high-frequency cut-offs of the Bk spectrum set the two break 
frequencies: the jitter frequency, ljj, being the main (low- 
frequency) spectral break and the high-frequency cut-off, cjf,, 
cf. Eqs. ( 15 I and (16 1. Below the jitter break, a; < Uj, 



the spectrum P{uj) is universal and flat. Between the breaks, 
ujj < to < (jJb, the spectral slope is solely uniquely deter- 
mined by the B^ spectral slope. Above tob, the spectrum is 
strongly suppressed. All these are in full agreement with the- 
oretical predictions, Eq. ( [T4| , [25|. Note that the radiation 
spectra are power-laws even though the particle distribution is 
monoenergetic. In contrast, the synchrotron spectrum exhibits 
an exponential fall-off above the peak frequency in this case. 
Note also that in isotropic turbulence, the low-frequency spec- 
trum is universal and flat, in contrast to the synchrotron rising 
spectrum with the power-law index being 1/3. 

Finally, we explored how the radiation spectrum is mod- 
ified when the radiation regime transits from small-angle to 
large-angle jitter. As 6 increases above unity, the spectrum 
attains some synchrotron features, namely, the peak is de- 
veloped and corresponds to the field-strength-dependent syn- 
chrotron frequency rather than loj and the low-energy slope 
tends to 1/3. We have demonstrated that different scales of 
magnetic turbulence contribute independently to radiation in 
different regimes: the large-scale part of B^ spectrum is in 
the large-angle regime, 5k ^ 1, whereas the small-scale part 
is still in the small-angle regime, 6k ^ 1. The resultant spec- 
trum is simply the sum of the two contributions. 

To conclude, the obtained results reveal strong inter- 
relation of transport and radiative properties of plasmas turbu- 
lent at sub-Larmor scales. They demonstrate how spectral in- 
formation can be a powerful tool to diagnose micro-turbulence 
in laboratory and astrophysical plasmas. 



In this paper we explored relativistic particle transport (dif- 
fusion) and radiation production in small-scale electromag- 
netic turbulence. We demonstrated that in the regime of small 
deflections, when the particle's deflection angle is smaller 
than the relativistic beaming angle, Aa ^ I/7, the pitch- 
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